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ABSTRACT 

We present an analysis of the binary and physical parameters of a unique pulsating 
white dwarf with a main-sequence companion, SDSS J1136-I-0409, observed for more 
than 77 d during the first pointing of the extended Kepler mission: K2 Campaign 1. 
Using new ground-based spectroscopy, we show that this post-common-envelope bi¬ 
nary has an orbital period of 6.89760103(60) hr, which is also seen in the photometry 
as a result of Doppler beaming and ellipsoidal variations of the secondary. We spectro¬ 
scopically refine the temperature of the white dwarf to 12 330±260K and its mass to 
0.601±0.036Mq. We detect seven independent pulsation modes in the light curve. 
A preliminary asteroseismic solution is in reasonable agreement with the spectroscopic 
atmospheric parameters. Three of the pulsation modes are clearly rotationally split 
multiplets, which we use to demonstrate that the white dwarf is not synchronously 
rotating with the orbital period but has a rotation period of 2.49 ± 0.53 hr. This is 
faster than any known isolated white dwarf, but slower than almost all white dwarfs 
measured in non-magnetic cataclysmic variables, the likely future state of this binary. 

Key words: asteroseismology, binaries: close, stars: white dwarfs, stars: individual 
(SDSS J113655.17-f040952.6) 


1 INTRODUCTION 

White dwarfs are the endpoint for all stars less massive than 
about 8 — 10 Mq , making them the fossils of more than 97% 
of Galactic stellar content (Williams et al. 2009; Doherty et 
al. 2015). They are representative of the general population 
of all stars, providing benchmarks that can calibrate theories 
of both stellar and binary evolution. 

White dwarfs do not sustain substantial nuclear burn¬ 


ing, and their evolution is almost exclusively defined by pas¬ 
sive cooling of the residual thermal energy of ions in the 
core. Fortunately, when hydrogen-atmosphere (DA) white 
dwarfs reach an appropriate effective temperature (roughly 
12, 500 K), their non-degenerate atmosphere develops a par¬ 
tial ionization zone, impeding energy transport and driving 
global pulsations; these are the ZZ Ceti stars (or so-called 
DAVs). Asteroseismology, enabled by matching the periods 
of observed luminosity variations to well-calibrated theoret¬ 
ical models, provides unparalleled insight into the interiors 
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Figure 1. Spectroscopic observations of the WD+dM binary J1136+0409. The top panel shows the discovery spectrum from SDSS. 
Underplotted in red is the best-fitting SDSS white dwarf model, as well as an M6 template in magenta estimating the contribution of the 
M-dwarf companion at redder wavelengths. The bottom left panel details our high signal-to-noise SOAR spectrum, from which we refine 
the white dwarf atmospheric parameters (red line, see Section 3). The bottom right panel shows two different epochs of our VLT/FORS2 
observations to calculate the radial velocity of the M-dwarf companion using the Nal doublet (see Section 4.1). 


of these white dwarfs (see reviews by Winget & Kepler 2008; 
Fontaine & Brassard 2008; and Althaus et al. 2010a). 

The development of a surface convection zone and the 
subsequent driving of pulsations are believed to be a natu¬ 
rally occurring phase for all white dwarfs (e.g., Bergeron et 
al. 2004; Castanheira et al. 2007). This so far holds regard¬ 
less of overall stellar mass (Giovannini et al. 1998) or core 
composition (Hermes et al. 2013b). 

White dwarfs found in binaries should also pulsate when 
they reach the appropriate temperature, and recently the 
first pulsating white dwarfs with main-sequence companions 
were discovered (Pyrzas et al. 2015). Roughly two-thirds of 
the known white dwarfs with M dwarf companions are wide 
enough that the progenitors evolved as if they were single 
stars (e.g., Nebot Gomez-Moran et al. 2011). The others 
evolved through a common envelope with a close companion 
and are post-common-envelope binaries (PGEBs). 

Most PGEBs are the detached precursors to cataclysmic 
variable systems. Nearly all progenitors of observed Super¬ 
novae la, independent of single vs. double-degenerate chan¬ 
nel, are believed to go through at least one common-envelope 
phase. Pulsating white dwarfs in PGEBs thus hold the ex¬ 
citing potential to allow for the first empirical test of the 
effects binary interaction has on the remnant white dwarf 
internal structure and chemical profiles, and can possibly 
directly constrain Supernovae la boundary conditions. 

Serendipitously, the only pulsating white dwarf in a 
confirmed PGEB discovered by Pyrzas et al. (2015), SDSS 
J113655.17+040952.6 (hereafter J1136+0409), was observ¬ 
able during the first science field of the two-wheel-controlled 


Kepler mission (K.S Gampaign 1). The extended Kepler mis¬ 
sion delivers space-based, time-series photometry on targets 
for up to 85 d along the ecliptic (see Howell et al. 2014). 

J1136+0409 was spectroscopically identified as a com¬ 
posite WD+dM from a Sloan Digital Sky Survey (SDSS) 
spectrum. Model-atmosphere fits to the Balmer absorption 
lines initially suggested the white dwarf had an effective tem¬ 
perature of 11 700 ± 150 K and log g = 7.99 ± 0.08 (Rebassa- 
Mansergas et al. 2012a), placing it in the empirical DAV in¬ 
stability strip. Pulsations were confirmed by a short (1.0 hr) 
run with the ULTRAGAM instrument mounted on the 3.5m 
New Technology Telescope, with significant periods near 
276.5 s and 182.2 s (Pyrzas et al. 2015). 

For nearly 78 d in 2014 August, the Kepler spacecraft 
observed J1136+0409 every minute, delivering an excep¬ 
tional light curve of the first known pulsating white dwarf 
in a PGEB. We present here a detailed analysis of the pul¬ 
sations and binary parameters of J1136+0409, enabled by 
space-based photometry, as well as spectroscopy from the 
8.2-m Very Large Telescope (VLT) and the 4.1-m Southern 
Astrophysical Research (SOAR) telescope. 


2 OBSERVATIONS AND REDUCTIONS 

2.1 Spectroscopic Observations 

The identification of J1136+0409 as both a candidate pul¬ 
sating white dwarf and a close WD+dM binary arose from 
an ancillary SDSS science fiber placed on the target. That 
spectrum, visible in the top panel of Fig. 1, clearly displays 
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Figure 2. A portion of the raw K2 Field 1 light curve obtained tor J1136+0409. This minute-cadence photometry has been extracted 
from the pixel-level K2 data and de-trended to correct for attitude tweaks in the telescope pointing, and has been 4cr clipped. The top 
panel shows 3d (less than 4 per cent of our full dataset) and is dominated by variability at the orbital period of 6.9 hr (see Section 4.2). 
The bottom panel shows a 4 hr stretch. The white dwarf pulsations are not visible above the point-to-point scatter. 1 ppt = 0.1 per cent. 


a composite of a white dwarf at bluer wavelengths and an M 
dwarf dominating at redder wavelengths. Fits for the white 
dwarf temperature put it within the empirical DAV insta¬ 
bility strip, and the SDSS sub-spectra showed significant 
radial-velocity variability over just 18 min, suggesting it is a 
close binary system. To better constrain both unique aspects 
of the system, we obtained follow-up spectroscopy. 

In order to solve the radial velocity of the dM compo¬ 
nent of the system, we obtained 11 service-mode observa¬ 
tions using the FORS2 instrument on the 8.2 m VLT/UTl 
(Antu) between 2009 April 4 and 2014 May 27. These 
intermediate-resolution spectra were all taken using the 
1028z grism and a 1" slit, covering a spectral region of 
7830 — 9570 A with a dispersion of roughly 0.8 A per pixel 
and 2.4 A resolution. We were exclusively interested in the 
sharp Nai absorption doublet at rest wavelengths of 8183.27 
and 8194.81 A (e.g., Rebassa-Mansergas et al. 2012b). Fur¬ 
ther details describing our reduction and analysis methods 
for the VLT data can be found in Schreiber et al. (2008). 

We display in the bottom right panel of Fig. 1 two 
epochs of the Nai doublet, separated by less than 2.0 hr 
on 2014 May 27. The analysis of these radial velocities and 
a full table of the observing epochs and exposure times are 
detailed in Section 4.1. 

Additionally, we sought to improve the white dwarf at¬ 
mospheric parameters, and obtained higher signal-to-noise 
(S/N) spectra using the Goodman spectrograph on the 4.1 m 
SOAR telescope (Clemens et al. 2004). These blue spectra, 
4x420s exposures taken consecutively on 2015 January 29, 
cover the higher-order Balmer lines and range from roughly 
3600 — 5300 A with a dispersion of roughly 0.8 A per pixel 
and 2.3 A resolution. We used a 3” slit, and the seeing was 
just under 1". Each exposure has a S/N-^ 70 per resolution 
element in the continuum at 4600 A. 

For all of our follow-up spectroscopy, we processed the 
images using the STARLINK packages FIGARO and KAPPA. 
The spectra were optimally extracted (Horne 1986) using 
the PAMELA package (Marsh 1989), and wavelength cali¬ 
brated using the many sky emission lines present. We sub¬ 
sequently flux calibrate and rebinned the spectra to a helio¬ 
centric frame using the molly package^. We have estimated 

^ http://www.warwick.ac.uk/go/trmarsh 


the spectral resolution from measurements of the full widths 
at half maximum of arc-lamp and night-sky emission lines. 

2.2 Space-Based Photometry 

The Kepler spacecraft provides an unprecedented capabil¬ 
ity for acquiring uninterrupted times series of faint variable 
stars, and the photometry of J1136-I-0409 (EPIC 201730811, 
Kp = 17.2 mag) opens a new era in seismology of white 
dwarfs in binary systems. Our short-cadence exposures, ob¬ 
served every 58.8 s, span nearly 78 days, from 2014 June 4 
03:02:00 UT to 2014 August 20 19:55:38 UT. 

With only two functional reaction wheels, the pointing 
of the Kepler spacecraft must be corrected periodically (with 
checks every 6 hr) by firing its thrusters. Thruster firings 
create significant discontinuities in the photometry and must 
be detrended. We extracted the light curve using a modified 
version of the kepsff code based on Vanderburg & Johnson 
(2014) within the pyke package maintained by the Kepler 
Guest Observer office. We have also applied an algorithm to 
select the pixel mask that looks at contiguous pixels above 
the background; our method is described in more detail in 
Hermes et al. (2014). A portion of the light curve is displayed 
in Fig. 2. 

Because the star has variability with a period near to 
6 hr, the thruster-firing period, we have utilized an inde¬ 
pendent extraction of the data using the tools outlined in 
Armstrong et al. (2015) to detrend the light curve. However, 
the r.m.s. scatter of this extraction was slightly higher, so 
we do not adopt it for our full analysis. 

We have removed all points with poor quality flags, as 
well as those affected by thruster firings, removing 32 265 
points. We subsequently clipped all points falling more than 
4x1.4826 times the median absolute deviation from the me¬ 
dian of the light curve, removing 351 points and ultimately 
leaving us with 91144 exposures over 77.8724 d. We have 
also excluded a brightening event 95.1 min in duration be¬ 
ginning at 2456844.36916 BJD(TDB); we have inspected 
the raw photometry using the k2flix package^ (Barentsen 
2015) and shown that this brightening corresponds to a solar 
system body passing through our target aperture. Pointing 

^ http://barentsen.github.io/k2flix/ 
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into the ecliptic, K2 encounters many asteroids as they pass 
through its field of view. 

Our final K2 light curve has a duty cycle of 79.7 per cent 
and has a formal frequency resolution of 0.149/rHz. The 
point-to-point scatter is evident by eye in Fig. 2 (expo¬ 
sures have a mean uncertainty of 13.0 ppt), and unfortu¬ 
nately overwhelms the ability to detect the pulsations by 
eye. Given the long baseline, the median level in a 1000- 
/iHz-wide, empty region of the Fourier transform centered 
at 1500 /iHz is roughly 0.064 ppt (64ppm). We describe am¬ 
plitudes in this manuscript in units of parts per thousand 
(ppt), where 1 ppt = 0.1 per cent relative amplitude. 


3 ATMOSPHERIC PARAMETERS 

From the original serendipitous SDSS spectrum of 
J1136-I-0409, Rebassa-Mansergas et al. (2012a) fitted the 
Balmer lines of H/3—He following the method of Rebassa- 
Mansergas et al. (2007). They found for the white dwarf a 
Teff = 11 700 ± 150 K and logg = 7.99 ± 0.08. 

We obtained higher-S/N SOAR spectroscopy to refine 
the atmospheric parameters of the pulsating component of 
J1136-I-0409. White dwarf atmospheric parameters are most 
sensitive to the higher-order Balmer lines, especially the 
surface gravity (e.g., Kepler et al. 2006), so we have used 
the blue-efficient Goodman spectrograph on the SOAR tele¬ 
scope (see Section 2.1). 

We have fitted the six Balmer lines H/1—H9 using the 
synthetic spectra computed by two independent groups. The 
first set of fits use the pure hydrogen atmosphere models 
and fitting procedure described in Gianninas et al. (2011) 
and references therein, which employ the ML2 /q; = 0.8 pre¬ 
scription of the mixing-length theory (Tremblay et al. 2010). 
We have fitted each of the four individual SOAR spectra, 
and found the weighted mean of the atmospheric parame¬ 
ters to be Teff = 12 680 ± 150 K and logg = 7.97 ± 0.03. 
In addition, we have individually fit the same spectra using 
the pure hydrogen atmosphere models detailed in Koester 
(2010), which also employ ML2/a = 0.8. From these fits we 
find a weighted mean of the atmospheric parameters to be 
Teff = 12 310 ± 150 K and log = 8.05 ± 0.03. 

To accommodate the unknown systematic uncertain¬ 
ties, we have chosen to adopt the average of these two in¬ 
dependent weighted mean fits, using the standard deviation 
of the two as the uncertainties: Tefi = 12 490 ± 260 K and 
log <; = 8.01 ±0.06. These values represent the effective tem¬ 
perature and surface gravity from fits to models with a one¬ 
dimensional treatment of convection. The SDSS fits used 
ML2/a = 0.6, which explains the disagreement with our 
new determination. 

Our adopted atmospheric parameters require a correc¬ 
tion to account for the higher-dimensional dependence of 
convection (Tremblay et al. 2013). Thus, our final best¬ 
fitting determinations are Teff = 12 330 ± 260 K and logg 
= 7.99 ± 0.06 which corresponds to an overall white dwarf 
mass of 0.601 ± 0.036 Mq and a radius of i?wD = 0.013 Rq 
using the mass-radius relations of Fontaine et al. (2001). 

By measuring the secondary contribution to the SDSS 
spectrum (see the top panel in Fig. 1), we estimate that 
the main-sequence component to J1136±0409 is of spec¬ 
tral type M6 ± IV. Using the empirical mass-radius relation 



Orbital Phase 


Figure 3. Radial velocity curve of the M dwarf companion 
in J1136-P0409, calculated from the Nal doublet. This curve is 
folded on the orbital period of 6.8976010 hr, and features a semi¬ 
amplitude of Asec = 222.6 ± 3.8 km s“^. 


Table 1. Journal of radial velocities from the Nal doublet. 


BJDxdb 


2452642.987095 

2452642.974837 

2452642.999318 

2454925.691559 

2454937.564039 

2456735.731791 

2456745.595020 

2456803.578535 

2456803.600854 

2456804.511249 

2456804.565577 

2456804.573015 

2456804.580532 

2456804.588222 


Radial Velocity 
(km s“^) 


76.6 ±56.1 
3.3 ±18.0 
124.9 ±26.5 
-118.7 ± 5.5 
-48.0 ±5.7 
-87.6 ± 11.0 
-78.0 ± 11.6 
-141.8 ±7.7 
-157.5 ±8.6 
-55.3 ±7.0 
192.3 ±6.1 

217.2 ±6.3 

235.2 ±6.4 
252.8 ±6.3 


Notes 


900s, SDSS subspec 
900s, SDSS subspec 
900s, SDSS subspec 
900s, VLT/FORS2 
900s, VLT/FORS2 
570s, VLT/FORS2 
570s, VLT/FORS2 
570s, VLT/FORS2 
570s, VLT/FORS2 
570s, VLT/FORS2 
570s, VLT/FORS2 
570s, VLT/FORS2 
570s, VLT/FORS2 
570s, VLT/FORS2 


of Rebassa-Mansergas et al. (2007), we estimate Msec = 
0.196 ± 0.085 Mq and Rsec = 0.195 ± 0.090 Rq. From the 
flux-scaling factors of each component, we can estimate the 
overall distance to the system (e.g., Rebassa-Mansergas et 
al. 2012a). From the white dwarf we find dwD = 125 ± 7 pc, 
and from the M dwarf we find dsec = 121 ± 56 pc, in good 
agreement. 


4 BINARY ANALYSIS 

4.1 Secondary Component Radial Velocities 

We have computed the radial velocities for each epoch of 
VLT/FORS2 and SDSS data using the same techniques de¬ 
scribed in Schreiber et al. (2008) and Rebassa-Mansergas et 
al. (2008). In summary, we fit a double-Gaussian-line pro¬ 
file to the Nai absorption doublet, with rest wavelengths of 
8183.27 and 8194.81 A, to compute a radial velocity for each 
epoch. The results are detailed in Table 1. We then compute 
a periodogram of the velocities (Scargle 1982). 

The folded radial-velocity curve of J1136±0409 is shown 
in Fig. 3. We have computed these parameters from a best 
fit to the equation 
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Orbital Phase 

Figure 4. The light curve of our K2 observations of J1136+0409 
folded at the RV-determined orbital period of 6.8976010 hr and 
repeated for clarity. We can produce a sensible model, overplotted 
in red, of Doppler beaming and ellipsoidal variations of the M 
dwarf that well describes the orbital photometric variability. 


where K^ec is the radial velocity semi-amplitude of the com¬ 
panion star, tfj is the time of inferior conjunction of the sec¬ 
ondary, Porb is the orbital period, and 7aec is the systemic 
velocity of the secondary. We find Porb =6.89760103(60) hr, 
7sec = 47.8 ± 3.5 km s“^, Tilsec = 222.6 ± 3.8 km s“^, and 
find for the inferior conjunction: 

fo(BJDTDB) = 2454925.5124(12) -f 0.287400043(25) E 

We have assumed circular orbits, and the best-fitting 
solution bears this sinusoidal signature. We find that the 
mass function is / = 0.328 ± 0.017 Mq. Adopting the spec¬ 
troscopic masses for both components, we find that the 
stars are separated by Ogep = 1-70 ± 0.08 Pq, and the 
volume-averaged secondary Roche lobe radius is roughly 
RLaec = 0.81 ± 0.07 Rq (Eggleton 1983). Thus, the M6 sec¬ 
ondary should sit comfortably within its Roche lobe. The 
system is far from the onset of mass transfer and its evo¬ 
lution into a semi-detached cataclysmic variable, which will 
occur in roughly 13 Gyr, when the system will have an or¬ 
bital period of 1.7hr (Schreiber & Gansicke 2003). 

4.2 Binary Light Curve Model 

The dominant feature of the light curve, visible by eye in 
Fig. 2, is the modulation near 6.9 hr. Taking a Fourier trans¬ 
form of the entire K2 dataset, we find a dominant peak at 
6.89756T0.00017 hr. We recognize this immediately as a less 
precise determination of the orbital period. 

We have phased the K2 photometry to the inferior con¬ 
junction defined from the radial velocities, and folded the 
light curve into 200 phase bins at the RV-determined orbital 
period. The result is displayed in Fig. 4. Both of our light 
curve extraction routines described in Section 2.2 produce 
an identical folded light curve. 

The asymmetric profile of the K2 light curve can be in¬ 
terpreted as a combination of reprocessed light (reflection), 
Doppler beaming, and ellipsoidal variations of the secondary 
component of J1136-1-0409. To confirm this, we have com¬ 
puted models using the lcurve software package developed 
by Tom Marsh, which computes models for binary systems 


containing at least one white dwarf (see Gopperwheat et al. 
2010). Our best-htting model is overplotted in red in Fig. 4. 

To calculate this model, we performed an initial explo¬ 
ration of the parameter space using a Monte Carlo anal¬ 
ysis, with initial guesses from the atmospheric parameters 
described in Section 3. This showed that the data is likely 
under-constrained, and that there are too many degeneracies 
to allow for convergence on a robust solution. We proceeded 
to fix essentially all parameters. 

We fixed the white dwarf radius and effective temper¬ 
ature as well as the secondary radius to the atmospheric 
parameters determined in Section 3, and used the corre¬ 
sponding Claret four-parameter limb-darkening coefficients 
of Gianninas et al. (2013) and Claret & Bloemen (2011). We 
convolved the best-fitting atmospheric model spectra (shown 
in the top panel of Fig. 1) with the Kepler bandpass to 
estimate the Doppler beaming factors for each component 
(e.g., Bloemen et al. 2011); we found — 2.34 and 

(^)aec ~ 11-69. The beaming factor for the M dwarf is quite 
high, and reflects the sharp spectral features that can move 
in and out of the Kepler bandpass. We hxed the mass ratio 
at q — Msec/MwD = 0.326, and assumed Kwo ~ qKsec- 

Our best fit has an inclination of i = 83 deg, for which 
the white dwarf contributes 61.6 per cent of the flux in the 
Kepler light curve. Aside from inclination, we also fit for the 
secondary temperature (which best fit as 3030 K) and the 
fraction of the irradiating flux absorbed by the companion 
(0.632). The fit is reasonable, with a Xred — 1-14. 

The slight bump around phase 0.25 in the light curve 
corresponds to ellipsoidal variations of the M dwarf, which 
peak twice per orbit with roughly the same amplitude; 
the second maximum near phase 0.75 is enhanced due to 
Doppler beaming. Ellipsoidal variations are tidal distortions 
of the secondary and are sensitive to the mass ratio, incli¬ 
nation, and secondary radius (e.g., Morris & Naftilan 1993). 
Unfortunately, our models are under-constrained without 
more information about the component masses (especially 
the primary radial velocity) or system inclination. From a 
lack of eclipses in the light curve, we can only constrain the 
inclination of the system to i < 85 deg. If we take the com¬ 
ponent masses from Section 3 at face value, the inclination 
must fall between 61 < i < 85 deg given the mass function. 

Doppler beaming imprints a modulation in the flux of 
a stellar binary as the source is approaching or receding, 
directly proportional to the radial-velocity semi-amplitude 
of the source as well as the slope of its spectrum through 
the observed filter (Loeb & Gaudi 2003; Zucker et al. 2007). 
The dominant source of Doppler beaming in J1136-1-0409 is 
the secondary, hence the peak nearest phase 0.75. The peak 
occurs slightly before phase 0.75 due to a small amount of 
reprocessed flux (reflection) from the much hotter primary. 

For now, we present our model as able to reproduce 
the data in hand, but do not yet consider it a robust, final 
solution. We have used our model to show that, with a fixed 
orbital period, the to from our raw, unfolded light curve 
is 2454925.5145(11) BJD(TDB), and agrees within Icr with 
the inferior conjunction defined from the spectroscopy in 
Section 4.1. 
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Figure 5. Fourier transforms of our K2 Campaign 1 observations of J1136+0409. The top panel shows the FT of the entire dataset out 
to the Nyquist frequency. The middle and bottom panels show the independent pulsation modes in more detail, marked with red dots; 
the full frequency solution is presented in Table 2. The dashed red line marks the significance threshold, described in the text. 


5 PULSATION ANALYSIS 

A Fourier transform (FT) of the 77.87 d K2 Campaign 1 
light curve of J1136+0409 is shown in Fig. 5. Table 2 details 
all of the significant periodicities detected in the light curve. 


5.1 Preliminary Mode Identifications 

The two significant peaks at lowest frequencies marked in 
blue in Fig. 5 fall at 40.2719 ± 0.0010 /rHz and its first har¬ 
monic of 80.5426 ± 0.0030 /rHz. These peaks correspond to 
the orbital frequency (/orb) and its first harmonic. Our light 
curve hts discussed in Section 4.2 reveal far more informa¬ 
tion about the binary parameters, so we will not discuss 
these signals further. 

Short-cadence Kepler data is contaminated by spurious 
frequencies arising from harmonics of the long-cadence ex¬ 
posure length (/lc = 566.479 ^Hz). We detect three such ar¬ 
tifacts with significance in our FT (7/lc, 8/lc, and 11/lc), 
and mark them in light gray in Fig. 5. We do not include 
these signals in our analysis. 

The pulsations of the white dwarf in J1136-1-0409 are 
confined to a region between 162.2 — 474.5 s, with a weighted 


Table 2. Significant frequencies in the K2 Field 1 photometry of 
SDSS J113655.17+040952.6 


ID 

Frequency 

Period 

Amplitude 

Phase 


(/rHz) 

(s) 

(ppt) 

(cycle) 

/i,+ 

3638.2374 

274.858366(88) 

3.493(50) 

0.1179(23) 

/l,0 

3578.5472 

279.44301(14) 

2.272(50) 

0.0547(35) 

/2 

5516.2411 

181.282868(73) 

1.841(50) 

0.4957(43) 

/s.o 

6164.0456 

162.231117(89) 

1.213(50) 

0.1201(65) 

/4,+ 

2961.1013 

337.71219(60) 

0.775(50) 

0.430(10) 

/4,- 

2848.1735 

351.10220(70) 

0.715(50) 

0.885(11) 

/5 

4955.8397 

201.78215(32) 

0.519(50) 

0.756(15) 

/3.+ 

6220.7698 

160.75181(24) 

0.448(50) 

0.641(18) 

/6 

2107.681 

474.4551(23) 

0.398(50) 

0.145(20) 

/3,- 

6107.362 

163.73682(27) 

0.399(50) 

0.515(20) 

fr 

2525.804 

395.9135(17) 

0.373(50) 

0.767(21) 

fh- 

3518.907 

284.1791(14) 

0.240(50) 

0.259(33) 

/orb 

40.2719 

24 831.21(62) 

4.063(50) 


2/orb 

80.5426 

12 415.79(47) 

1.347(50) 


8/lc 

4531.8157 

220.66210(20) 

0.994(50) 


7/lc 

3965.3057 

252.18737(56) 

0.467(50) 


11/lc 

7930.6581 

126.09294(18) 

0.354(50) 
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mean period of 263.3s (Clemens 1993). Such relatively 
short pulsation periods are entirely consistent with the 
spectroscopically determined temperature of 12 330 K (see 
Mukadam et al. 2006); hot DAVs near the blue edge of the 
instability strip have shallower convection zones and thus 
drive shorter-period pulsations than do their cooler coun¬ 
terparts (Goldreich & Wu 1999). 

We have computed a significance threshold following 
the method outlined in Greiss et al. (2014). In short, we 
have randomly shuffled the fluxes for each point in the light 
curve, and computed the highest peak in the entire resultant 
FT for 10 000 random permutations. The 2a and 3cr lines 
are defined by the lines for which 95.5 and 99.7 percent 
of the highest overall amplitude recorded in each shuffled 
light curves fall below that threshold. They are 0.330 and 
0.359 ppt, respectively. The threshold we adopt is conserva¬ 
tive since the basis light curve includes significant signals, 
notably the orbital modulation. However, no questionable 
signals fall between our 2cr and 3a lines, so we adopt the 
99.7 percent threshold (0.359 ppt) for simplicity. 

We identify 16 signals in the FT above the significance 
threshold, which we describe in Table 2. The quoted period, 
amplitude, and phase uncertainties in Table 2 come from a 
simultaneous least-squares fit to the data using the software 
package Period04 (Lenz & Breger 2005). 

All signals cleanly prewhiten from the light curve. The 
only frequency with a marginally significant residual occurs 
at 2961.585^tHz (with 0.341 ±0.050 ppt amplitude). We have 
also prewhitened the light curve by three very low frequen¬ 
cies, at 0.11842 ^Hz, 0.6624 ^Hz, and 0.8118/iHz, respec¬ 
tively. We do not believe these signals are astrophysical in 
nature, but are rather long-term trends that have not been 
perfectly removed in the light curve extraction. 

We note that since the white dwarf is contributing 
roughly 62 per cent of the flux in the Kepler bandpass, the 
pulsation amplitudes are all underestimated by a factor of 
roughly 1.61. The amplitudes have not been adjusted in Ta¬ 
ble 2, but this should be taken into account when comparing 
to pulsation amplitudes observed at other wavelengths. 

In order to perform an asteroseismic analysis, we must 
know the spherical degree {£) of the pulsation mode in ques¬ 
tion. We begin with the initial assumption that all modes are 
the two lowest degrees possible {i = 1,2) out of simplicity 
and because £ ^ 3 modes suffer from significant geometric 
cancellation (Dziembowski 1977). Additionally, we note that 
three regions of the FT have signals separated in frequency 
by a similar amount, roughly 57/rHz. We assume that these 
are rotationally split multiplets (see Dolez & Vauclair 1981), 
and discuss these signals further in Section 5.3. 

In Table 2 we have identified seven independent pul¬ 
sation modes, three of which we identify as having more 
than one rotationally split component. We have tentatively 
identified the central (0), prograde (±), and retrograde (—) 
components of these three modes. In one case, fi,-, we have 
used the expectation that this may be a rotationally split 
component of a triplet to relax the significance threshold, 
and include this low-amplitude mode in our analysis. 

All of the pulsation signals we have identified in Table 2 
are extremely stable in both amplitude and phase. Such fre¬ 
quency stability is a feature of hotter DAVs, which show 
extremely stable pulsation modes (e.g., G117-B15A, Kepler 
et al. 2005). In fact, every mode in J1136±0409 shows excep- 


Table 3. Summary of the observed periods and the best-fitting 
asteroseismic model. We have fixed the I marked by a * symbol. 


ID 

Observed 
^obs (®) 

Model 

^calc (®) 

£ 

k 

Ck,i 

/i 

279.44 

280.78 

1* 

3 

0.350 

/2 

181.28 

183.41 

2 

4 

0.054 

/3 

162.23 

164.77 

1* 

1 

0.490 

/4 

344.28 

345.15 

1* 

5 

0.493 

/5 

201.78 

199.73 

2 

5 

0.157 

/e 

474.46 

471.19 

1 

7 

0.412 

fr 

395.91 

396.23 

2 

12 

0.113 


tional amplitude and phase stability, such that if we break 
the K2 data into four ~3-week subsets, a linear least-squares 
fit for all values are consistent within the statistical uncer¬ 
tainties. 


5.2 Preliminary Asteroseismic Analysis 


We have used the seven independent pulsation modes iden¬ 
tified from the K2 Gampaign 1 light curve as input for 
an initial asteroseismic analysis. Our asteroseismic analy¬ 
sis proceeded in the same way as described in more detail 
in Bischoff-Kim et al. (2008). 

In short, we matched the observed pulsation periods 
(fobs) to theoretical periods (Pcaic) computed from adiabat- 
ically pulsating a grid of models derived from an updated 
version of the White Dwarf Evolution Gode (Lamb & van 
Horn 1975). Our core composition profiles follow the smooth 
models of Safaris et al. (1997), and feature a sharp transi¬ 
tion between the carbon and helium layer, where diffusive 
equilibrium is assumed for this transition and no diffusion 
coeffecients are used. 

The models had five variable parameters: effective tem¬ 
perature, overall white dwarf mass, the mass of the hydro¬ 
gen envelope (Mh), the central oxygen abundance relative 
to carbon (Vo), and the location of the edge of the homo¬ 
geneous carbon-oxygen core (Afm). Our model grids were 
spaced in steps of 200 K, 0.005 Mq, 0.10 dex, 0.10 dex, and 
0.05 dex, respectively, with a range of 10 000 — 13 000 K, 
0.55 - 0.80 Mq, 10"'‘“ - 10"® °, 0.50-0.99dex, and 0.1- 
0.8 dex, for a total of 1541202 models. The mass of the 
helium envelope (Mhb) was fixed at MHe/Af*= 10~^ ®. 

We assume the rotationally split modes with three com¬ 
ponents correspond to a pulsational triplet, and thus iden¬ 
tify /i, /a, and /4 as ^ = 1 modes. We make no assumption 
about the identification of the other modes, and minimized 
the overall residual ctrms: 


ctrms = 


E(-Pcalc - Pobs)" 

'^obs 


( 2 ) 


where riobs = 7 is the number of observed modes. 

The best-fitting asteroseismic solution is detailed in Ta¬ 
ble 3, with a CTRMS = 2.02 s. It arose from the following 
parameters: Tefj^ 12000 K, Mwd = 0.605 M©, Mh/M* = 
10“^®, Xo = 0.99, and Xtm = 0.70. The top 20 best fits 
all shared the same £ identifications, with 2.09 < ctrms < 
2.47 s: their means and standard deviations form an aster¬ 
oseismic temperature (12260 ± 270 K) and mass estimate 
(0.605 ± 0.014 Mq). The other parameters were consistent 
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Figure 6. Contours showing the goodness of fit, (jrmS: from our asteroseismic analysis of J1136+0409. The scale runs linearly in tenths 
of seconds; darker regions correspond to a better fit. The best asteroseismic solution is marked as a white dot and the la spectroscopically 
determined atmospheric parameters are marked in blue. This is a five-dimensional parameter space, so to make these two-dimensional 
projections the best-fitting values for the other parameters were optimized at each grid point. 


within 0.1 dex of the best-fitting asteroseismic solution. 
The shaded contours in Fig. 6 show only solutions with 
CRMS < 3.5 s, and represent just 0.2 per cent of all models. 

The best-fitting models consistently prefer a high 
oxygen-carbon ratio in the core, but expectations of this 
parameter are fairly unconstrained (see the discussion in 
Bischoff-Kim & Metcalfe 2011), and Xo has a relatively 
small impact on DAV period spectra (see Table 5 in Bischoff- 
Kim et al. 2008). The models of Althaus et al. 2010b predict 
Xo ~ 0.7 for a 0.6 Mq white dwarf. 

Given how precisely we were able to determine the pul¬ 
sations periods with our K2 light curve (to an average of 
0.00066 s), the best-fitting asteroseismic model has a rather 
large crms of 2.02 s. This indicates that our best fit may not 
have fully converged on a final solution, or that the models 
inherently lack some physical details. However, we are en¬ 
couraged by the good agreement between the atmospheric 
parameters (temperature and mass) determined by both 
asteroseismology and spectroscopy. Our preliminary aster¬ 
oseismic analysis also improves our estimate of the white 
dwarf rotation period. 


5.3 White Dwarf Rotation Period 

The symmetric frequency spacing of the /i, /a, and /r 
modes, illustrated in Fig. 7, suggests that these are rotation- 
ally split multiplets of three independent pulsations modes. 
Since we see no evidence that these modes have five com¬ 
ponents, we assume that these are rotationally split triplets 
of three f = 1 multiplets, each of a different radial overtone 
number {k). 

To first order, the frequency splitting is related to the 
rotation of the white dwarf by the relation 

Siy = m(l — Ck,e)^ (3) 

where Q is the rotation frequency of the star and the Ck,e 
value represents the moment of inertia of the mode (Unno 


et al. 1989). In the asymptotic limit, Ck,e ~ l/f’(^ + !)• 
However, this does not hold for trapped modes, or modes of 
low radial order {k < 10), which is very likely the case for the 
relatively hot DAV in J1136-I-0409. Fortunately, our best¬ 
fitting asteroseismic model has predictions for these Ck,e 
values, which we have listed in Table 3. Notably, /i has 
a very strong deviation from the asymptotic value (0.5). 

We find the frequency splitting for each of the three 
triplets: <5/i = 59.6890(29) /rHz, Sfs = 56.7057(92) /rHz, and 
5/4 = 56.464(11) ^Hz. For each triplet, we computed the 
spacing between the central component and the prograde or 
retrograde component, and we average these values weighted 
by their uncertainties. For /4 we have simply taken half the 
frequency difference between / 4 ,_ and f4,,+ - Using the Ck,i of 
our best asteroseismic model, the observed frequency split¬ 
tings correspond to rotation periods of 3.02 hr, 2.50 hr, and 
2.49 hr. 

Given the excellent agreement for /s and / 4 , we do not 
use the trapped mode /i in our estimate of the overall rota¬ 
tion period but rather use it to constrain the uncertainties 
in our determination. We conclude that this white dwarf is 
rotating at 2.49 ± 0.53 hr. 

One troubling aspect of the three rotationally split 
triplets is the inconsistent relative amplitudes: the central 
component of /s has the highest amplitude of the multiplet, 
whereas the central component of /4 is not detected. The 
central component of fi has the second-highest amplitude, 
while the prograde component appears strongest. If the sys¬ 
tem is at high inclination, as we found in Section 4.2, geo¬ 
metric considerations predict the amplitude of m = ±1 are 
much larger than m = 0 for i = \ modes, and that m = ±2, 0 
are favoured for i = 2 modes (see Pesnell 1985). The ex¬ 
pected amplitude ratios work for fi, but not the other two 
multiplets. We note that if the multiplets are all m = ±2, 0 
components of £ = 2 modes rather than m = ±1,0 compo¬ 
nents of £ = 1 modes, the splittings would correspond to a 
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Figure 7. Detailed regions of the Fourier transform of 
J1136+0409 showing the three different rotationally split pulsa¬ 
tion modes, all viewed with the same frequency scale. The triplets 
are split by roughly 57 /rHz, which corresponds to a rotation pe¬ 
riod of roughly 2.49 hr (see Section 5.3). Our interpretation sug¬ 
gests that the central (m = 0) component of /4 is not observed. 


rotation period of 8.3 ± 1.8 hr using the asymptotic relation 
for f = 1, 2 modes. 

Stochastically excited stars, such as solar-like oscilla¬ 
tors, have equipartition of energy, such that the amplitudes 
of all components of a rotationally split multiplet are sym¬ 
metric and predictable (Gizon & Solanki 2003). However, 
observationally, this does not appear to be the case for pul¬ 
sating white dwarfs: different multiplets in the same star 
almost never have the expected relative amplitudes for a 
given inclination and the same intrinsic mode amplitudes. 
For example, the amplitudes of both the 1=1 and the 
1 = 2 modes of one of the best-studied pulsating white 
dwarfs, PG1159—035, show inconsistent relative amplitudes 
(Winget et al. 1991). Similarly, the amplitude ratios of an¬ 
other class of pulsating compact objects, hot subdwarfs, do 
not show consistent relative amplitudes for all = 1 modes 
of the same star (Pablo et al. 2012). Thus, the different rela¬ 
tive amplitudes of the observed rotationally split multiplets 
in J1136+0409 is not problematic to our analysis. 

We have searched the FT for a signal at the correspond¬ 


ing spin frequency, /spin = 111.6 /rHz, but do not detect any 
peaks in this region above 0.2 ppt. We do not see a photo¬ 
metric signal of the white dwarf rotation period. 

Finally, any slight asymmetry in what should be the 
even frequency splitting of a pulsation multiplet can con¬ 
strain the presence of a weak magnetic field, which also acts 
to break spherical symmetry and cause subtle changes in the 
frequency spacing (Jones et al. 1989). We find that the two 
triplets with all three components, fi and /s, are symmetric 
within the uncertainties except for the prograde component 
of /i, which is still symmetric within 2cr. Thus, we can only 
put an upper limit on the magnetic held strength. Follow¬ 
ing the method and the scaling factor detailed in Winget 
et al. (1994), the frequency splitting symmetry to better 
than 0.023 /rHz roughly constrains the dipolar magnetic held 
of the white dwarf in J1136+0409 to < 10 kG. Such a low 
global magnetic held may be consistent with the lack of a 
corresponding photometric signal of the white dwarf spin. 


6 DISCUSSION 

The pulsating component in J1136+0409 is the hrst carbon- 
oxygen core white dwarf known to pulsate and also have 
a close, detached companion (Pyrzas et al. 2015). The 
only other pulsating white dwarfs with close companions 
(which have also evolved through at least one and probably 
two common-envelope phases) are the extremely low-mass 
(Mwd < 0.25 Mq) white dwarfs showing pulsations (Her¬ 
mes et al. 2012, 2013a,b; Kilic et al. 2015), as well as the no- 
longer-detached pulsating white dwarfs in cataclysmic vari¬ 
ables (e.g., van Zyl et al. 2004). 

As discussed in Section 5.2, the seven independent pul¬ 
sation modes in J1136+0409 are exceptionally well deter¬ 
mined, and provide the possibility for the first in-depth anal¬ 
ysis of the interior structural effects of common-envelope 
evolution. Unfortunately, our asteroseismic analysis is ham¬ 
pered by the fact that there are nearly as many free param¬ 
eters in our models as there are observed pulsation modes. 
We have shown that the 20 best-fitting asteroseismic models 
are in good agreement with the spectroscopically determined 
atmospheric parameters. However, the best fit is not excel¬ 
lent: it has a ctrms = 2.02 s, which is significantly higher 
than the uncertainty for each pulsation mode (cr < 0.001 s). 

For these preliminary asteroseismic fits, we have con¬ 
sistently found a hydrogen layer mass, which controls the 
rate at which the white dwarf cools, that is thinner than 
we expect for an isolated white dwarf. Our best fit finds 
this layer has a mass of 10“'*'®, thinner than the 

canonical value for DA white dwarfs, Mh/M*= 10“'^ ° (e.g., 
Althaus et al. 2005). The seismologically derived hydrogen 
layer mass is nearly an order of magnitude smaller than ex¬ 
pected from shell burning in the evolutionary models, and 
could indicate that some of the hydrogen mass was stripped 
in the common envelope phase. However, the fact that we 
have not considered thinner helium-layer masses and fixed 
that parameter at the maximal value of 10“^'° 

may have an influence on the derived hydrogen-layer mass. 

The only previous constraints on the hydrogen-layer 
mass of a PCEB have required comparing precision stud¬ 
ies of deeply eclipsing systems to model-dependent mass- 
radius relationships. Such studies have shown that the white 
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dwarfs in PCEBs contain both a relatively thick hydrogen- 
layer mass, Mh/A/*« 10“^, as in the case of the 0.535 Mq 
NN Ser (Parsons et al. 2010), as well as a relatively thin 
hydrogen-layer mass, Mh/M*« 10 “®, as in the case of the 
0.439M 0 SDSS J1212-0123 (Parsons et al. 2012). 

Pulsations also have the potential to further refine the 
binary parameters of this unique system. As mentioned in 
Section 5.1, the pulsations visible in J1136-I-0409 are notably 
stable, and a single frequency can well-describe each peak in 
the FT. However, we know that the phases of the pulsations 
are periodically changing at the orbital period. This is a 
light-travel-time effect as the white dwarf moves around its 
centre of mass (e.g., Barlow et al. 2011). 

Constraining this phase variation would allow us to di¬ 
rectly measure Kwd and provide a model-independent so¬ 
lution to the physical parameters. Given the mass estimates 
from Section 3, we expect the amplitude of this phase varia¬ 
tion tober = 1.8±0.7s. However, the K2 light curve is sim¬ 
ply too noisy to enable our directly measuring this quantity. 
We also attempted, unsuccessfully, to directly detect side- 
lobes in the amplitude spectrum, manifest as splittings from 
the central components of our pulsations set apart by the 
orbital frequency. We expect the sidelobes to have roughly 
2 per cent the amplitude of the central component, which is 
simply buried in the noise of the FT (Shibahashi & Kurtz 
2012). We can only set an upper limit of r < 10 s on this 
phase variation from the K2 data. 

Our radial-velocity observations were carried out to ex¬ 
clusively cover the Nai absorption doublet at wavelengths 
longer than 7500 A, so we have no spectroscopic informa¬ 
tion on the white dwarf radial velocities. Measuring the 
white dwarf radial-velocity semi-amplitude directly from 
the Balmer lines should be possible, since in all likelihood 
Kwd > 70 km s“^. This will fully solve the orbital param¬ 
eters, yielding a precise estimate for the mass ratio, as well 
as a gravitational redshift for the white dwarf. 

Despite our inability to use the oscillation phases to 
constrain the binary mass ratio, the exceptional stability of 
the pulsations can be monitored longer-term to watch the 
orbital period evolution of the system. Roughly 90 per cent 
of detached, eclipsing PCFBs observed for more than 5yr 
show significant changes in the mid-eclipse times incon¬ 
sistent with secular evolution (for a list, see Zorotovic & 
Schreiber 2013). Some of these eclipse variations have been 
attributed to light-travel-time effects caused by circumbi- 
nary planets (e.g.. Marsh et al. 2014). However, follow-up 
studies have been unable to directly detect these external 
companions and have been damaging to the circumbinary in¬ 
terpretation of eclipse timing variations (Hardy et al. 2015). 

With follow-up photometry it will be possible to use 
the stable pulsations in J1136-1-0409 as a clock to moni¬ 
tor for external light-travel-time effects. Modulations in the 
pulse arrival times must come from light-travel-time effects, 
whereas mid-eclipse-time variations can arise from orbital 
period variations caused by other means, such as Apple¬ 
gate’s mechanism (Applegate 1992). To enable further mon¬ 
itoring, we have included in Table 2 the phase point for the 
time-of-maximum of each pulsation mode from the K2 ob¬ 
servations. These phases align with the mid-point of the K2 
light curve, 2456851.5625920 BJD(TDB). 

Our K2 observations of the white dwarf in J1136-1-0409 
provide one of the most precise constraints on the rotation 


period of an object that has gone through common-envelope 
evolution and is still fully detached. The only other well- 
constrained rotation period of a PCEB has been determined 
for the eclipsing system V471 Tau, which has a 9.25-min 
spin period, although the K2V companion is nearly Roche- 
lobe filling and may have previously transferred mass and 
angular momentum onto the white dwarf (see Hussain et al. 
2006 and O’Brien et al. 2001). 

In the context of other asteroseismically deduced ro¬ 
tation periods of white dwarfs, J1136-I-0409 is rotating 
more rapidly than any other known isolated white dwarf 
with clearly resolved pulsation multiplets (Kawaler 2004, 
2014). The previous record-holder is the hot pre-white-dwarf 
PG 2131-1-066, which rotates with a period of 5.07 hr; no¬ 
tably, it is also in a wide binary with a K7V companion, 
although the two are currently separated by at least 200 an 
and so did not undergo a common-envelope event (Kawaler 
et al. 1995; Reed et al. 2000). 

Over time, gravitational radiation and magnetic brak¬ 
ing from the M dwarf will bring the two components of 
J1136-I-0409 closer together, and the system will likely be¬ 
come a cataclysmic variable. Interestingly, all non-magnetic 
white dwarfs in cataclysmic variables have rotation periods 
on the order of minutes or faster, much more rapid than ob¬ 
served for the white dwarf in J1136+0409 (e.g., Sion et al. 
1994, 1995; Szkody et al. 2002). Over time, mass transfer in 
cataclysmic variables is expected to transport enough angu¬ 
lar momentum to spin up the white dwarf to near break-up 
velocities (e.g., King et al. 1991). 

For now, the components of J1136+0409 are detached 
and still well separated, and the secondary fills roughly 24 + 
20 per cent of its Roche lobe. It is therefore very unlikely that 
episodes of stable mass transfer occurred in the past that 
could have considerably affected the current white dwarf 
rotation. Therefore, the measured rotation period (2.49 ± 
0.53 hr) and the upper limit on the magnetic field (< 10 kG) 
put interesting constraints on angular momentum transport 
in evolved stars. Notably, angular momentum transport in 
giants is still an open problem in astrophysics (Cantiello et 
al. 2014). First constraints on the problem of core spin-up 
and envelope coupling in red giants are being explored using 
seismic data on first-ascent giants from the original Kepler 
mission (Beck et al. 2012; Mosser et al. 2012). 

If isolated white dwarf rotation periods reflect some sort 
of coupling with the red-giant envelope, then the relatively 
fast rotation of the white dwarf in J1136+0409 could have 
resulted from truncated red-giant evolution via common- 
envelope ejection. Additionally, the relatively low magnetic 
field strength of the white dwarf in J1136+0409 may have 
precluded the development of magnetic torques to slow the 
white dwarf rotation after the common-envelope event (Suijs 
et al. 2008). 

Finally, we have shown that the white dwarf rotation 
period is not synchronized with the orbital period of roughly 
6.9 hr. Similar results from three hot subdwarfs with M 
dwarf companions in the original Kepler mission showed a 
similar lack of synchronization, although those subdwarfs ro¬ 
tate much slower than the orbital period (Pablo et al. 2011, 
2012). Finding more such systems help constrain theoretical 
models for timescales of tidal synchronization. 


© 0000 RAS, MNRAS 000, 000-000 


Internal effects of common-envelope evolution 11 


Table 4. Summary of system parameters determined for 
J1136+0409. Parameters were derived from (1) VLT/FORS2 
spectroscopy, (2) SOAR/Goodman spectroscopy, (3) SDSS spec¬ 
troscopy, and (4) K2 Campaign 1 light curve. 


Parameter 

Value 

Source 

Orbital period, Porb 

6.89760103(60) hr 

(1) 

WD temperature (spectroscopy) 

12 330(260) K 

(2) 

WD temperature (seismology) 

12 260(270) K 

(4) 

WD logg 

7.99(06) cgs 

(2) 

WD mass (spectroscopy) 

0.601(13) Mq 

(2) 

WD mass (seismology) 

0.605(14) Mq 

(4) 

WD radius (seismology) 

0.013 Rq 

(4) 

Secondary spectral type 

M6 ± IV 

(3) 

Secondary radius 

0.195(90) Rq 

(3) 

Secondary mass 

0.196(85) Mq 

(3) 

WD distance 

125(07) pc 

(3) 

Secondary distance 

121(56) pc 

(3) 

Secondary RV, Kaec 

222.6(3.8) km s”! 

(1) 

Orbital separation, asep 

1.70(08) Rq 

(1,2,3) 

Secondary roche lobe 

0.81(07) Rq 

(1.2,3) 

WD rotation period 

2.49(53) hr 

(4) 

WD magnetic field 

< lOkG 

(4) 

WD hydrogen layer mass 

Mh/M*^^ lO-'^-® 

(4) 


7 SUMMARY AND CONCLUSIONS 

We have used extensive ground-based spectroscopy and 
space-based photometry to constrain the physical parame¬ 
ters of the first known pulsating white dwarf that underwent 
a single common-envelope event, J1136+0409. We summa¬ 
rize these values in Table 4. 

Both our time-series spectroscopy and photometry vary 
at 6.89760103 hr, which we determine is the orbital period 
of this evolved binary. Our photometry from the extended 
Kepler mission shows orbital variability that can be well- 
reproduced with a model of Doppler beaming, reflection, and 
ellipsoidal variations of the MOV main-sequence companion. 
We have matched the seven observed independent pulsation 
modes to theoretical evolution models, and find reasonable 
agreement between the spectroscopically determined white 
dwarf temperature and mass. Asteroseismology also allows 
us to peer into the interior of the white dwarf and constrain 
its hydrogen layer mass to 10“^ ®, although we 

only consider this a preliminary asteroseismic investigation 
and urge caution in over-interpreting this result. 

Notably, the observed splittings of the pulsations 
demonstrate that the white dwarf in this close binary is 
rotating at a period of 2.49 ± 0.53 hr, more rapidly than the 
orbital period. These splittings are symmetric within the 
uncertainties, which sets an upper limit on the dipole mag¬ 
netic field of the white dwarf to < 10 kG. Both values are 
extremely useful inputs for further modeling of the previ¬ 
ous common-envelope event that shaped the present binary 
configuration of J1136+0409, as well as angular momentum 
transport in stars that have likely underwent truncated red- 
giant evolution. 
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